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We study a system of two coupled kicked rotors, both classically and quantum mechanically, 
for a wide range of coupling parameters. This was motivated by two published reports, one of 
which reported quantum localization, while the other reported diffusion. The classical systems are 
chaotic, and exhibit normal diffusive behavior. In the quantum systems, we found different regimes, 
depending on the strength of the coupling. For weak coupling, we found quantum localization similar 
to that exhibited by single kicked rotors. For strong coupling, we found a quasi-diffusive growth of 
the width of the momentum distribution, in which the apparent diffusion coefficient decreased as time 
increased. The behavior in this regime can be described by the scaling theory of weak localization 
for two-dimensional disordered lattices. The weak and strong coupling regimes are separated by 
a regime of complex intermediate behavior. Thus we explain the apparent contradiction in the 
literature as being due to qualitatively different regimes of behavior, which we call strong and weak 
quantum localization. 

PACS numbers: 05.45.Mt, 05.45.Pq, 72.15.Rn, 05.45.Ra 



I. INTRODUCTION 

Kicked rotor models [l|, 0, Q have played a prominent 
role in the study of classically chaotic systems and their 
corresponding quantum behavior. Although quite sim- 
ple, these models exhibit most of the distinctive features 
of classical and quantum chaos. The particular phe- 
nomenon that is the subject of this paper is quantum 
localization, which is a localization in angular momen- 
tum space that is in sharp contrast to the diffusive motion 
typical of classical chaos. The understanding of quantum 
localization is considerably facilitiated by the existence of 
a mapping from the rotor model to a model of a parti- 
cle moving in a disordered lattice Q , which exhibits the 
phenomenon of Anderson localization 0]. 

Single rotor systems have been thoroughly studied and 
are well understood However, the situation is less 
clear for systems of two coupled kicked rotors. Indeed, 
there are examples in the literature whose results appear 
to be contradictory. Doron and Fishman Q have studied 
one such model, and they found quantum localization of 
the Anderson type, scaling exponentially with the cou- 
pling strength. On the other hand, Adachi et al. 0] have 
studied a model with different coupling potentials, and 
have obtained results that exhibit diffusive growth in the 
width of the state. They concluded that the coupling 
between the two rotors can restore the mixing that 
would otherwise be suppressed by quantum localization. 

There are several possible explanations for this appar- 
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ent contradiction. The two papers are based on mod- 
els that differ in the form of the kicking and coupling 
potentials. Although it is logically possible for the two 
models to behave qualitatively differently, this would be 
very surprizing, since the forms of the potentials are sim- 
ilar enough that we should expect them to belong to the 
same generic class. The ranges of interaction parameter 
strength used in the two papers are quite different. In- 
deed, the lattice emerging from the choice of parameters 
in [rj turns out to be periodic, which would give rise to ex- 
tended Bloch states. So we need to study the two models 
in equivalent ranges of paramenters. The two papers also 
used different criteria for identifying localization. Doron 
and Fishman look for an exponential fall-off in the 
angular momentum distribution, while Adachi et al. Q 
look for a saturation in time of the width (standard de- 
viation) of the distribution. So we need to apply both 
criteria to both models to determine whether they agree. 
A further complication arises because two-dimensional 
disordered lattices contain the marginal case of weakly 
localized states, which creates more interesting possibili- 
ties. 

In this paper we investigate these questions. To do so, 
we consider a moregeneral Hamiltonian that includes the 
models of @ and Q as special cases. We study the sys- 
tem numerically for a suitable range of parameters, and 
examine the shape and fall-off of the angular momen- 
tum distribution, together with its standard deviation 
as a measure of the state width. In Section II we dis- 
cuss the Hamiltonian, the initial state, and the method 
of numerical calculation. In Section III the classical be- 
havior of the system is studied, and is shown to lie within 
the chaotic region of phase space for the chosen range of 
parameters. Section IV presents the numerical results for 



2 



the quantum system. In Section V we examine whether 
the results are consistent with the scahng properties of 
weak locahzation. Section VI discusses the conclusions. 



momentum representations is achieved by Fast Fourier 
Transform (FFT) and its inverse. The same process is 
then repeated for the next period, and so on. 
The initial state is chosen to have the form 



II. MODEL 

The general Hamiltonian is chosen to be 
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T{pi,P2) = "1 Pi + ^ "2 P2 



V {61,62) = [Ai cos 61 + A2 cos 62 + A3 cos 61 cos 61 

+XiCos{ei-62)] (1) 

The dimensionless parameters are related to the physical 
quantities as follows: 
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Here r is the time span between kicks, and /i, I2 are 
the moments of inertia of the first and second rotors, 
respectively. The parameters Ai and A2 are the single- 
rotor kicking strengths, while A3 and A4 are the couplings 
between the two rotors. 

The choice Ai = A2 = A4 = 0.0 and ai = 1, 012 = 
\/2 , h~l will yield the Hamiltonian studied in Q . The 
main systems analyzed in Q correspond to A3 = 0.0, 
«! = a2 = ^ X 27r, and ai = 02 = x 2tt. Since 
these values of a are commensurate with tt, they yield 
periodic lattices, giving rise to extended Bloch-like states 
and excluding the possibility of Anderson localization. 
Therefore, we have used a different value in our numerical 
studies. 

The quantum system is obtained by replacing the 
momentum variables with operators pi = —id/d6i and 
P2 = —id 1 362- The time evolution for one period of the 
kicking Hamiltonian is given by the Floquet operator 
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The numerical calculation is performed in two stages. 
Starting in the momentum representation, the propaga- 
tor e~*'^(Pi'P2\ which is diagonal in this basis, is applied 
for the duration between two kicks. Then the kicking 
propagator ^-i-^i^^^^'i) is applied in the angle representa- 
tion. The transformation between the angle and angular 
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TABLE L Parameters of the initial state 
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Rotor no. 


7710 


Am 


So 


1 





1.25786 





2 





1.25786 






This is a product of Gaussian wave packets, centered 
at momentum toq ^md angle 6q, with momentum width 
Am = 1.25786. Hence we have SAm = ?i/2. A sum- 
mary of the parameter values are given in Table [D 

The momentum values form a discrete grid. The angle 
values are also set on a grid of the same size, since 
the Fast Fourier Transforms are instances of the discrete 
Fourier transform algorithm. The grid size is chosen to be 
2" = 2048 grid points, ranging from m = -1023 to 1024 
of dimensionless angular momentum values. The numer- 
ical calculation is performed for 30,000 kicking steps. The 
coupling parameters A3 and A4 range from 0.0 to 3.0 
in steps of 0.5, for each value of Ai and A2 in the set 
{0.0, 0.25, 0.5}. The simulations are done for the values 
oi\ ^ 0L2 = 1.0, which are not commensurate with tt, 
thus ensuring that the corresponding Anderson lattice is 
non-periodic. 



III. CLASSICAL DIFFUSION 

In order to analyze the quantum dynamics, it is neces- 
sary to determine whether the classical counterpart lies 
in the chaotic regions of phase space. Only when it is ver- 
ified that the corresponding classical evolution is chaotic, 
can the localization in the quantum system be positively 
attributed to quantum effects. It is also interesting to 
compare the behavior of the quantum system to its clas- 
sical counterpart. 

The classical mapping of each point (6'["' , 6*2"' , Pi"' , P2"' ) 
of the phase space from its values at kick n to kick n + 1 
is: 
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For the classical calculations, an ensemble of 1,000,000 
initial states was created, with the same Gaussian distri- 
bution in phase space as the initial quantum state. Each 
member of the ensemble was evolved separately using the 
mapping ([5]) . The results of the calculation show similar 
chaotic behavior for all of our choices of parameters, with 
the exception of Ai = A2 = A3 = and A4 > 0, for which 
the total angular momentum is a constant of motion. The 
variance of the angular-momentum distribution for each 
rotor exhibited diffusive growth with time. 



{Pfl - {P2f = D2 t 
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where the average is over the classical ensemble. This 
is illustrated in FiglTJ From this we conclude that the 
chosen parameters do lie within the classically chaotic 
region of parameter space. 
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FIG. 1: (Color online) Width of the angular momentum dis- 
tribution for a classical rotor [defined similar to ®] a: Ai — 
0.5, A2 = 0.5, A3 = 3.0, A4 = 3.0. b: Ai = 0.5, A2 = 0.5, A3 = 
1.0, A4 = 2.0. c: Ai = 0.5, A2 = 0.5, A3 = 1.0, A4 = 0.0. 



For our chosen initial states, we have (p) = in all 
cases, so we can replace the variance with (p^)- The de- 
pendence of the diffusion coefficients, Di, on the kicking 
and coupling parameters can be analyzed in the lowest 
order (so-called quasi-linear) approximation by assum- 
ing a uniform probability distribution for the angle after 
a sufficiently large number of kicking steps n, and inte- 
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This yields 
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A comparison between this approximation and the 
numerically computed diffusion constant is shown in 
FiglH The computed diffusion coefficients oscillate above 
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FIG. 2: Ratio of the classical difltusion coef. to the quasi-linear 
approx., Di/Di, vs couplings A3 and A4, for Ai — 0.5, A2 = 
0.5. 

and below the approximate value. This is similar to a 
known result for the standard map, for which one has 
= where K is the kicking strength in the stan- 
dard map. This result was also verified by a more sys- 
tematic method similar to that used in Q. 



IV. RESULTS FOR THE QUANTUM ROTORS 

The evolution of the quantum system yields a 
momentum-space wave function, ^'(^'1,^2)7 ^s a func- 
tion of time, for each set of the Hamiltonian parameters. 
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From it, we obtain the two-rotor momentum probability 
distribution, |^(pi,P2)Pi and calculate the widths (stan- 
dard deviation) of the single-rotor momentum distribu- 
tions, 



S^^{{p~{p)f) (9) 

as a function of time (measured in kick numbers). These 
arc plotted as a function of the square root of time, rather 
than time, in Figures [31 HI and because this choice 
of variable is effective in showing both short and long 
time behavior on the same scale. This plot would yield 
a straight line for normal diffusion. 
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FIG. 3: (Color online) Localized regime, 
a: Ai = 0.5, Aa = 0.5, A3 = 2.5, A4 = 0.5. 
b: Ai = 0.25, Aa = 0.25, A3 = 2.5, A4 = 0.5. 
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FIG. 5: (Color online) Quasi-diffusive regime, a: Ai — 
0.5, Aa = 0.5, A3 = 3.0, A4 = 3.0. b: Ai = 0.25, Aa = 
0.25, A3 = 3.0, A4 = 3.0. 



Quasi-difFusive: The width of the momentum distribu- 
tion increases with time, but the growth rate di- 
minishes, as illustrated in FiglH 

Intermediate: The width of the momentum distribu- 
tion may exhibit a complex oscillatory behavior, as 
illustrated in FiglH 

Tables |lT] and IIIII show the occurence of the three be- 
haviors as a function of the coupling strengths A3 and A4, 
for fixed values of the kicking strengths Ai and A2. For 
vanishing kicking strengths, Ai = A2 = 0, all states with 
non-zero coupling strengths appeared to be localized. 
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FIG. 4: (Color online) Intermediate regime, a: Ai — 0.5, Aa — 
0.5, A3 = 2.5, A4 = 1.0. b: Ai = 0.25, Aa = 0.25, A3 = 
2.5, A4 = 1.0. 

We observed three kinds of behavior of the momentum 
width as a function of time. 

Localized: The width of the momentum distribution 
saturates in time, as illustrated in Fig[3) 



TABLE II: Behavior of the angular momentum distribution 
width for fixed kicking parameters Ai = Aa = 0.50. D: Quasi- 
diffusive regime. I: Intermediate regime. L: Localized regime. 



A4 


0.0 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


A3 

0.0 


L 


L 


L 


L 


L 


L 


L 


0.5 


L 


L 


L 


L 


L 


I 


I 


1.0 


L 


L 


L 


L 


I 


I 


I 


1.5 


L 


L 


L 


I 


I 


I 


D 


2.0 


L 


L 


I 


I 


I 


D 


D 


2.5 


L 


L 


I 


I 


D 


D 


D 


3.0 


L 


I 


I 


D 


D 


D 


D 



The quasi-diffusive behavior occurs for the largest val- 
ues of the couplings between the rotors. Figure [5] shows 
that at the earliest times, the quantum rotor behaves dif- 
fusively, as does the classical rotor, but at later times the 
quantum diffusion rate decreases. It is not clear whether 
at very long times the quantum rotor will continue to be- 
have diffusively but with a reduced diffusion coefficient, 
or whether it will eventually saturate at some large local- 
ization scale. We shall return to this question in the next 
section. 
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TABLE III: Behavior of the angular momentum distribution 
width for fixed kicking parameters Ai = A2 = 0.25. D: Quasi- 
diffusive regime. I: Intermediate regime. L: Localized regime. 



A4 


0.0 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


A3 

0.0 


L 


L 


L 


L 


L 


L 


L 


0.5 


L 


L 


L 


L 


L 


L 


L 


1.0 


L 


L 


L 


L 


L 


L 


L 


1.5 


L 


L 


L 


L 


L 


L 


L 


2.0 


L 


L 


L 


L 


I 


I 


I 


2.5 


L 


L 


L 


I 


I 


I 


D 


3.0 


L 


L 


L 


I 


I 


D 


D 
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FIG. 7: (Color online) Localized regime: angular momentum 
g distribution for rotor #1 along the section p2 = 0. Ai = 




tion at early times, a: classical rotor, b: quantum rotor. -600 -400 -200 200 400 600 

Ai = 0.5, A2 = 0.5, A3 = 3.0, A4 = 3.0. angular momentum 



The two localization criteria - exponential fall-off of 
the angular-momentum distribution, and saturation of 
the momentum-distribution widths - are consistent with 
each other. For states that are classified as "Localized" 
in the tables, we have verified that after sufficient time 
has elapsed, the momentum distribution falls off expo- 
nentially with a rate that is independent of time. For 
the "Quasi-diffusive" states, the momentum distribution 
continues to spread. This is illustrated in Figures [7] and 
[SI repectively, where we plot the angular momentum dis- 
tribution for pi along the section P2 = 0. 

V. SCALING THEORY AND THE 
QUASI-DIFFUSIVE REGIME 

In the Quasi-diffusive regime, the computations do not 
clearly determine whether the motion will remain diffu- 
sive at very long times, or whether it will become local- 
ized at some scale that is greater than the grid size of the 
computation. In this section we apply the scaling theory 
of Anderson localization (see [l^l for a review) to answer 
this question. 

In its usual context, the phenomenon of Anderson 



FIG. 8: (Color online) Quasi-diffusive regime: angular 
momentum distribution for rotor #1 along the section P2 ~ 0- 
Ai = 0.5, A2 = 0.5, A3 = 3.0, A4 = 3.0. a: at t=1000, b: at 
t=5000, c: at t=30000 steps. 



localization refers to a charged particle moving in a dis- 
ordered medium. Scaling theory considers the behavior 
of the conductance g (not the conductivity) as a func- 
tion of the system size L. The mean- free-path I may 
be defined as the distance beyond which the phase of the 
wave function is essentially randomized by the scattering. 
No phase coherence exists between regions separated by 
more than I, and scaling behavior is possible only for L 
greater than I. The conductance go — g{l) a.t this cutoff 
length acts as a measure of the disorder in the system. 

It is useful to distinguish the regimes of strong and 
weak disorder. Strong disorder gives rise to exponen- 
tially localized states, and to a non-ohmic conductance 
that scales as g{L) oc exp(— L/^), where C is a localization 
length that is generally greater than I, but may approach 
I in extreme cases. In d = 1 dimension all states exhibit 
this strong localization, while in d = 3 dimensions this 
strong localization occurs only for sufficiently strong dis- 
order. The case of d = 2 dimensions may exhibit strong 
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localization, but it can also exhibit a marginal localization 
for weak disorder due to the coupling of the two rotors, 
for which g{L) differs from ohmic conduction by a weakly 
varying correction proportional to log(L). Weak disorder 
usually leads to ohmic conductance, however d = 2 is an 
exceptional case. 

Abrahams et al [ll[ have argued that the logarithmic 
derivative of the conductance g with respect to sample 
size L is a function I3{g) of the conductance alone, 



d\ogg 
dlogL 



Pig) 



(10) 



The form of P{g) for ohmic conduction in d dimensions 
is 



Pi9) = {d-2) 
For strong localization, it takes the form 
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Pig) = log 



go 



(11) 



(12) 



There exists a critical value for go at which the large scale 
behavior changes from conducting to localized. 

For weak disorder, it is possible to use perturbation 
theory to calculate the corrections to ohmic conduction, 
which yields 



Pig) ^ id -2) 



g 



(13) 



where a is some constant for the system. It is clear that 
d = 2 is a special case. 

We can now increase the size of the system L (L > I). 
From (fT3ll for d = 2, we have 



dlogL 



a 

g 



(14) 



Integrating this with respect to L from the lower bound 
I yields 



giL) = ga-a log 



(15) 



The localization length, (, is defined as the length scale 
L for which g{L) = 0, that is 



(16) 



C = I exp 

In transport theory for weak scattering [l^l, we have 

go oc I (17) 

so we can rewrite ()16p as 

C = le^^ (18) 
where & is a constant for the system. 



Recall that the kicked rotor system can be mapped 
onto an Anderson- type lattice [j]. If the phase of the 
rotor's wave function is randomized between two kicks, 
then the average displacement in momentum after the 
first kick can be regarded as the counterpart of the mean- 
free-path in the Anderson lattice. This should be true 
for the large values of the rotor coupling parameters, for 
which we obtained quasi-diffusive behavior. The average 
displacement in momentum after the first kick is 



irl + 



rl) 



o,r\ 



ri ,r2 — — oo 



where 



f/o,r (0,0|f/|ri,r2), 



iri,r2) 



(19) 



(20) 



and U is the Floquet operator defined in ([3]). The matrix 
elements are in the angular-momentum basis. If, in the 
quasi-diffusive regime, the state is actually localized on 
some large scale, then we may expect the scaling theory 
of Anderson localization to apply. 

The initial state for our numerical calculations is a 
Gaussian, effectively localized in a finite portion of the 
infinite angular-momentum lattice. The support of the 
wave function grows with time. At any finite time, the ef- 
fect of lattice sites far outside the wave function's support 
will be negligible, and it would not make any difference 
to the dynamics if the wave function were instead located 
on a finite lattice of a size not less than the width of the 
state. So it is reasonable to regard the state width S (or 
some multiple of it) as being the analog of the sample 
size L in the Anderson lattice. 

It is also reasonable to assume that the diffusion in 
angular-momentum space of the rotor will contain the 
same kind of information as does the diffusion of charged 
particles on the Anderson lattice, and that they will scale 
similarly with S or L, respectively. According to the Ein- 
stein relation [l^ , the diffusion coefficient is proportional 
to the mobility, and hence to the conductivity. But in 
d = 2 dimensions, the conductance g scales the same 
way with size as does the conductivity. Thus the diffu- 
sion coefficient D of the rotor should scale the same as 
does the conductance g of the lattice. 

The measure of the angular-momentum state width of 



the rotor at time t is taken to be S = 



ip-{p)y 



By analogy with the classical equation ^ , we define the 
diffusion coefficient D of the kicked rotor to be the slope 
of the line relating to t. If the relationship is not a 
straight line, as is the case for Fig.® and ^ (which, 
however show S vs Vi), then we can define a time- 
dependent diffusion coefficient as the local slope of the 
curve of S'^(t) vs t. 

In the spirit of the scaling hypothesis, we assume that 
the time-dependent diffusion coefficient D at any time 
i is a function only of S at that time, that is to say. 
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D = D{S(t)). In analogy with JTH), we write 
dlogD 



d logs' 



a 
D 



(21) 



with a being some unknown constant for the system. 

In the kicked-rotor model, all quantities change only 
at the discrete time steps t„. The state width at time tn 
is Sn- From the definition of Z)„ = D{tn) as the slope of 
the curve of S'^ vs t, we have 



£.2 



S„ 



D„AT 



(22) 



where AT — t„+i — i„. According to the scaling assump- 
tion, we should have a relation of the form 



Dn = cl — a log 



T 



(23) 



by analogy with (jlSp . with two as yet undetermined para- 
meters c and a. Hence we obtain 



AT cl - aloe 



Sn 

T 



(24) 



As in the original scaling argument that introduced the 
localization length in the Anderson lattice, we can de- 
fine the saturation length A as the value of Sn that sets 
the right hand side of the above equation to zero. This 
yields a result analogous to (fT5|) . 



(25) 



with b = c/a. Since the corresponding parameter b was 
a constant for the Anderson lattice, we expect that here 
b will be approximately constant in the Quasi-diffusive 
regime of the rotor if the condition of marginal localiz- 
ation for weak disorder is valid. 

We now wish to test the scaling theory by comparison 
with the results of numerical simulation. In the previous 
section, we computed S^ as a function of t, for a range of 
the kicking and coupling parameters. These results are 
now regarded as data. It should be noted that, although 
the diffision coefhcient D = d{S'^)/dt plays a fundamen- 
tal role in developing the scaling theory, it is not neces- 
sary to compute that derivative numerically. Only the 
data S'^ vs t is needed for the computation. 

To carry out the scaling-theory calculation, we first 
compute / from ([TO]) . The matrix elements in can 
be obtained by running the numerical simulation for one 
kicking step from an initial state that is localized at pi ~ 
P2 — 0. These values of / are contained in Tables HVl and 
IVl Next we use the recurrence relation to obtain 
the scaling theory values for Sf-^, which we denote as S*^, 
using Sq = P as the starting value. The parameters c 
and a cannot be determined directly, so we must do the 
recursion calculation for several values of c and a, and 
determine the best values of these two parameters by 
minimizing the sum of squared differences between the 
scaling-theory values S^ and the previously computed 
data for 5^(t„). 



TABLE IV: Numerical values of /, b and A [See (O and ((25]l ] 
for each run with Ai = A2 = 0.50 







A4 


1.0 


1.5 


2.0 


2.5 


3.0 


A:, 
1.0 










3.57 


1.5 








3.37 


3.86 


2.0 






3.20 


3.67 


4.15 


2.5 




3.06 


3.52 


3.98 


4.46 


3.0 


2.96 


3.39 


3.84 


4.30 


4.77 



E 



A4 


1.0 


1.5 


2.0 


2.5 


3.0 


A3 

1.0 










0.87 


1.5 








0.86 


1.01 


2.0 






0.84 


1.02 


1.19 


2.5 




0.97 


1.00 


1.11 


1.11 


3.0 


0.82 


0.93 


1.14 


1.22 


1.16 



A4 


1.0 


1.5 


2.0 


2.5 


3.0 


A3 

1.0 










80 


1.5 








61 


190 


2.0 






47 


155 


579 


2.5 




60 


115 


330 


630 


3.0 


34 


79 


306 


816 


1207 



The choice of AT in and must be made so 
that S'^ (t) varies nearly linearly with t during the interval 
AT. Taking AT = 1 kicking step has the disadvantage of 
yielding a curve whose slope might fluctuate considerably 
from step to step. It is better to choose a larger value 
of AT so as to smooth out such short-time fluctuations, 
but of course AT must not be so large as to cover a 
range over which the curve is appreciably nonlinear. We 
found the choice of AT = 300 kicking steps to be a good 
compromise, but the results are insensitive (in the first 
decimal place) to the precise value that was chosen. 

The least-squares fitted values of a and c were used 
to calculate the parameter b = c/a. The values for b 
and the corresponding saturation length A, from (|25p . 
are listed in Tables IIVI and [V] The values of b seem to 
become nearly constant as we move to the lower right 
of the tables, away from the transition zone between the 
Localized and Quasi-diffisive regions shown in Tables |ll] 
and mil This is what we should expect if the scaling 
hypothesis is valid, and our Quasi-diffusive regime is act- 
ually a regime of marginal localization for weak disorder. 
In contrast to the nearly constant parameter b, the local- 
ization length A varies over orders of magnitude in these 
ranges of kicking and coupling strengths. Our results are 
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TABLE V: Numerical values of I, b and A for each run with 
Ai = A2 = 0.25 







A4 


1.0 


1.5 


2.0 


2.5 


3.0 


A3 

1.0 










3.54 


1.5 








3.34 


3.83 


2.0 






3.17 


3.64 


4.13 


2.5 




3.03 


3.49 


3.96 


4.43 


3.0 


2.92 


3.36 


3.81 


4.27 


4.75 



E 



A4 


1.0 


1.5 


2.0 


2.5 


3.0 


A3 

1.0 










0.60 


1.5 








0.65 


0.71 


2.0 






0.57 


0.67 


0.80 


2.5 




0.57 


0.69 


0.79 


0.86 


3.0 


0.58 


0.67 


0.77 


0.83 


0.79 







A4 


1.0 


1.5 


2.0 


2.5 


3.0 


A3 

1.0 










30 


1.5 








29 


58 


2.0 






19 


42 


112 


2.5 




17 


39 


90 


201 


3.0 


16 


32 


72 


149 


202 



consistent with marginal localization for weak disorder 
on the scale of A. 



VI. CONCLUSIONS 

This paper was motivated by an apparent contra- 
diction between the conclusions of two published pa- 
pers [1,01 on very similar models of coupled kicked rotors. 
One group found quantum localization, while the other 
reported diffusive behavior. It was unclear whether the 
differing conclusions were due to the different criteria for 
localization used by the two groups, or whether they re- 
flected fundamental differences between the two models, 
or some other reason. By studying a more general model, 
of which those two are special cases, we were able to re- 
solve the apparent contradiction. The two models do, 
indeed, belong to the same generic class and have simi- 
lar behavior. When the two criteria for localization are 
applied to the same model in the same range of para- 
meters, they always agree. The differences between the 
results reported by the two groups were primarily due 
to their use of very different interaction strengths and 
kicking intervals. 



Beyond resolving this discrepancy in the literature, our 
study of two coupled kicked rotors with two different cou- 
pling interactions has revealed some interesting system- 
atic behavior as the kicking and coupling strength para- 
meters are varied. The classical model exhibits diffusive 
motion. An analytic formula for the dependence of the 
classical diffusion coefficient on the kicking and coupling 
strength parameters was derived, by means of a random- 
phase approximation. The numerically computed diffu- 
sion coefficient oscillates about this approximate formula, 
much as occurs for uncoupled rotors. As the couplings 
between the rotors become large, the behavior of the ro- 
tors is determined mainly by the sum of the two coupling 
parameters, A3 + A4. This generalization holds, at least 
qualitiatively, for the quantum rotors too. 

The quantum mechanical model shows classical diffu- 
sion at short times, but at long times the diffusion is lim- 
ited by a kind of quantum localization. For zero coupling 
strength, the model reduces to two independent kicked 
rotors, each of which exhibits the well-established quan- 
tum localization that is the analog of Anderson localiz- 
ation in a one-dimensional random lattice. This single 
rotor localization behavior persists for a moderate range 
of couplings between the two rotors. 

As the coupling strength increases further, the system 
passes through a complex transition zone, before entering 
a quasi-diffusive regime for strong coupling. In this quasi- 
diffusive regime, the mean squared width of the momen- 
tum distribution increases with time, but the slope of 
the curve (the local diffusion coefficient) gradually de- 
clines. Within practical limits of grid size and simulation 
time, it was not possible to determine directly whether 
the motion is bounded on a much larger scale or grows 
without bound. We found that the marginal localization 
for weak disorder due to the coupling of the two rotors 
fits well to this quasi-diffusive motion, and thereby we 
conclude that the state is very probably localized on a 
larger scale than our practical computational grid. We 
were also able to estimate this large localization scale, 
which varies smoothly with the coupling strength para- 
meters. In the intermediate region between strong and 
weak localization, the state-width oscillates in a compli- 
cated way, but its envelope always lies beneath that found 
in the regime of weak localization. 

Thus we have shown that a system of two coupled 
kicked rotors exhibits a complex set of behaviors, includ- 
ing a transition from strong to weak localization. This 
would be expected from the formal mapping of the two- 
rotor system onto a two-dimensional Anderson lattice. 
However, this correspondence is not trivial, since the 
theory of Anderson localization assumes a random lat- 
tice potential, whereas the effective potential generated 
by the kicked-rotor system is not random, but merely 
quasi-periodic with incommensurate frequencies. Thus 
the existence of strong and weak localization in the rotor 
system is not guaranteed by the theory of localization in 
random lattices. 

However, we disagree with one conclusion of Adachi et 
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al 0], who claimed that, for sufRciently large coupling 
strength, the coupled quantum rotors exhibit normal dif- 
fusion. Their results were obtained for a special set of 
parameters, for which the dimensionlcss kicking interval 
([2]) is exactly commensurate with tt. Such special cases 
give rise to periodic lattices, for which the quantum states 
arc not localized, regardless of the kicking and coupling 
strengths. These highly special cases do not correspond 
to typical physical situations. It is more realistic to study 
the generic non-commensurate case. We believe that, in 
the generic case, the two coupled quantum rotors exhibit 
weak localization on a long time scale. 

Finally, we mention a possible connection between this 
model and experiment. Experimental realization of sys- 



tems like the one considered in the present paper are 
nowadays intesnly explored. The experimental observa- 
tion of 2-dimensional optical localization have been re- 
ported [l^, [iBl , where the localized results were shown 
to match with the predictions of scaling theory. 
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